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We study the long wavelength limit of a spin 1/2 Heisenberg antiferromagnetic two-leg ladder, 
treating the interchain coupling in a non-perturbative way. We perform a mean field analysis and 
then include exactly the fluctuations. This allows for a discussion of the phase diagram of the system 
and provides an effective field theory for the low energy excitations. The coset fermionic Lagrangian 
obtained corresponds to a perturbed SU(4)i/U(l) Conformal Field Theory (CFT). This effective 
theory is naturally embedded in a SU (2)2 x Z2 CFT, where perturbations are easily identified in 
terms of conformal operators in the two sectors. Crossed and zig-zag ladders are also discussed 
using the same approach. 

I. INTRODUCTION 

Antiferromagnetic Heisenberg spin ladders have been a subject of central interest during the last years. These 
are intermediate systems between the gapless critical spin 1/2 Heisenberg chain and the ordered spin 1/2 2D system 
relevant for undoped cuprate superconductors. The simplest realization, i.e. the two-leg ladder, shows a dramatic 
different excitation spectrum with respect to the one of an isolated chain. It has a finite gap to the first excitation and 
magnetic correlations are short ranged. Several inorganic compounds have been recently synthesized and modeled as 
Heisenberg ladders . Exponential decay of the low temperature magnetic susceptibility was the first signal of the 
existence of a spin gap in two-leg ladder materials. Neutron and optical measurements also manifest the presence of a 
gap and are consistently described by a two-leg ladder model with exchange integrals of the same order in the chains 
direction(J) and along the rungs (J'). 

Theoretically, the existence of a gap was early predicted from numerical exact diagonalization and strong coupling 
perturbation theory (J/ J' <C 1) More recently field theoretical technics have been used to analyze the excitation 
spectrum in the weak coupling regime (J'/ J <C 1) ||f§. These treatments give access to the whole low energy 
excitation spectrum as well as to the dynamical susceptibilities, which are essential to compare with experimental 
probes. The philosophy underlying this study is the following: spin operators are expressed in the well known bosonized 
representation of each chain and the interchain coupling is treated as a small perturbation in this representation. The 
applicability of these studies is then valid in principle only in the weak coupling regime and its use in the description 
of e.g. the experimentally realized two-leg ladders in which J' ~ J/2 should be taken with some care. It is therefore 
not clear up to which value of the results of |3|,|| are applicable, and it is important to develop theoretical methods 
which could be used beyond the weak coupling regime. 

The picture that emerges from the weak coupling analysis leads to a description in terms of triplet of massive 
Majorana fermions and a singlet Majorana fermion with a different mass (which has been estimated to be minus three 
times the triplet mass) 0. The only interactions between these fermions are marginal current-current terms which 
have been argued to simply renormalize their masses and velocities. A question, that has been risen in recent studies 
of the Raman scattering spectrum [^J, is whether marginal interactions can in fact be disregarded. In particular, 
correlation functions obtained disregarding marginal interactions apparently do not fit experiments (see e.g. ||). 

In this work we analyze the complete phase diagram of the two-leg antiferromagnetic ladder. Our approach, first 
used here for spin ladders, starts from a fermionic representation of the spin operators in the functional integral 
framework, as introduced in for spin chains. With a simple Ansatz to the Mean Field (MF) configurations we 
show that the system undergoes a cross-over from a weak to a strong coupling regime at an intermediate value of 
J'/ J. We then introduce fluctuations around MF and take them into account to all orders to construct the low energy 
effective field theory. 

The resulting theory corresponds to a coset Conformal Field Theory (CFT) of symmetry SU(A)i/U(l)i SO , perturbed 
by relevant operators (of dimension 1) and marginal operators (of dimension 2) arising from the single occupancy 
constraint as well as from the amplitude fluctuations of the link fields introduced to decouple the fermionic interactions. 
It should be noted that our approach is based on the assumption that the local single occupancy constraint can be 
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implemented as a very last step, while it is taken into account globally from the beginning. The correctness of this 
procedure is not guaranteed from first principles, but is supported a posteriori. 

We show that the complete structure of these perturbations can be retained and that they take a simple form in 
the language of conformal embeddings. In particular, the marginal terms which arise can be easily classified in the 
new language and their effect can then be studied in a non-perturbative way. When the relevant perturbations are 
expressed in the embedded SU (2)2 x Z2 CFT language the spectrum is naturally separated in the triplet and singlet 
of Majorana fermions. These results, which are valid up to J'/ J « 8/ir 2 , extend to finite coupling the weak coupling 
study of ||]. It should be stressed that recent estimates of the ratio of exchange constants lead to values of J' / J 
around 1/2 in several cuprate materials |l0| . 

In order to illustrate the generality and ease of use of our approach, it is then applied to the so-called crossed 
ladders and zig-zag ladders. Phase diagrams and low energy theories are obtained in the region containing the weak 
coupling limit; further analysis and details will be considered elsewhere. 

The paper is organized as follows. In Section II we introduce the model, present Hubbard-Stratonovich decoupling 
technics and perform a MF analysis, discussing the resulting phase diagram. In Section III we construct the low 
energy effective field theory: our theory contains four Dirac fermionic species corresponding to the spin and band 
indices of the ladder. In Section IV we show that the theory has a natural relation to SU {2)2 x Z2 CFT through 
conformal embedding (the last part arises from the two electronic bands). In Section V we briefly report results on 
crossed and zig-zag ladders. Finally, in Section VI the conclusions and possible further developments of our method 
are given. 



II. MEAN FIELD ANALYSIS 



We consider the Heisenberg Hamiltonian for a two-leg spin 1/2 ladder, 

N 2 / V \ 

H = EE Wft + yS« +1) (1) 

n=l 1=1 V Z J 

where N is the number of sites along the chains, J > and J' > are the couplings between adjacent spins along the 
legs and rungs respectively. For mathematical convenience we assume periodic boundary conditions (P.B.C.) in both 
directions (notice that the Hamiltonian is suitable written for arbitrary n-leg ladders; in the present case the physical 
coupling along the rungs is in effect J'). 

The spin variables can be represented in terms of fermionic operators with spin Cn,a as 

where a are Pauli matrices, together with a local constraint that ensures one spin per site, imposed on the physical 
states by 

ct«lc$ a \phys) = \phys) . (3) 

Throughout this paper we will not use the summation convention neither for site nor leg indices; repeated spin (Greek) 
indices are summed. 

Using (||) and (0) the Hamiltonian ([j]) can be rewritten as 

h = - ( -c^c {l) c m c {l) + — c t(o c a+i) c t(i+i) c (o \ +c (4) 

n=l 1=1 ^ ' 

where C — —N/2 is an irrelevant constant term. 

We now trade (|4|) for a quadratic Hamiltonian via a Hubbard-Stratonovich transformation, at the usual price of 

introducing auxiliary fields Bn associated to terms containing c^f\ a c^ a and B' n associated to terms containing 

c n l p c n8- ft ^ s na tural to interpret as localized on the leg (I) links between sites n and n + 1 and B' n as localized 
on the rung links. After the transformation the Hamiltonian reads 
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n=l 1=1 
,, N 2 

+ t E E (^k^ 15 + ^iTMi + «) (5) 



As we look for a low energy effective theory, we treat the B variables in a long wave approximation. To this end, 
we parameterize these fields in terms of real MF values (Bo, B ) and fluctuations 

B« = B exp(iaA® + aS®), B n = b' exp(iaA n + aR n ) . (6) 

Notice that we have included both phase and amplitude fluctuations, which will play important different roles in the 
following. For this reason, we explicitly distinguish the Hermitean (Rn , Rn ), and anti-Hermitean (iA$ , iA^p) parts 
of the fluctuation fields. The expression for B n will be eventually modified when B — (see eq. (|3l|)). 

As a first step, we perform the MF evaluation of the Hamiltonian (0) by setting the fluctuations to zero. The 
resulting MF Hamiltonian is then a tight-binding model for two coupled chains, 

JV 2 N 

„t( 2 U! 



H , - -t W (c^c (l) + c t(0 c (0 "l - 2t' V (JWcW 

n— 1 1=1 n— 1 



C n,/3 C n,0 



4iV 2 8iV ,2 

where 

, J B . J'B 

— > * — ■ (8) 

The coupled tight-binding model is easily diagonalized by means of a double Fourier transform. We first decouple 
two bands by means of 

c (1) = — fc (+) -c (_) ) (9) 

v/2 n,a!i \ u ) 

c (2) = — fr (+) +c ( - ) ) (10) 
and then introduce pseudo-momentum operators d^ + \d^' by 



r (+) 



1 w o 

— £d« exp(-i — ), (11) 

m— 1 



1 w o 

^E^exph— ), (12) 



in terms of which the Hamiltonian reads 

AT 



H mf = -J2 (2i cos(^m) + 2t')4«d« 



N 

m— 1 



N 

^ (2*cos(^m) - <L + ^ 2 + ^i' 2 - (13) 



m— 1 



This expression clearly represents a decoupled two-band tight-binding model. 

The constraint (|J), meaning one electron per site, forces the system to be exactly at half filling. Low-energy 
excitations are then achieved by creating holes just below the Fermi surface and creating electrons just above it |TT| . 
Notice that this can be done only if 
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It'l < 1*1, 



(14) 



that is when the Fermi level crosses both bands. If this condition is not satisfied, the system presents a finite energy 
gap to spin excitations. 

The actual values of t and t' are determined by minimizing the energy of (|l3|). In order to perform this evaluation we 
introduce a lattice spacing a and a position coordinate x = na (x G [0, L = Na\); the appropriate pseudo- momentum 
coordinate is k = m2ir/{Na) (fc £ [— 7r/a, 7r/a]). The mean-field Hamiltonian then reads 

H mf = -^- f ' {2tcos(ka) + 2t')d^+\k)d^\k)dk- 



2tt 



■k / a 



l r /a 



2tt 



— n/a 



4i 



(2icos(fca) - 2t')d^- ] (k)d { ~ ] {k)dk 
8L + , 2 



-t ^ 1 

aJ aJ 

from which we can read the dispersion relations for each band, sketched in Fig. 1, 

e (+) (fc) = -2tcos(ka) - 2t' , 



(15) 



(16) 



£ (-)(fc) = -2tcos(fca) + 2t'. 



(17) 




FIG. 1. Dispersion relation for the two-band tight-binding model 

The Fermi momentum for each band is defined through the equation 

e(+)(4 + ))= e (-)(4-)). (18) 

The local constraint in eq. ^) leads to the global constraint N 1 -^ — 2N, where is the occupation number 

operator for each band. Besides, = ^^k^ . We thus obtain 



(19) 



these implying 



4 +) = */a-k£\ 



5(+)(4+)) = e (-)(4 _) ) = o. 



(20) 



(21) 



The values of t and t' are now determined by minimizing the energy of (llq) under half-filling conditions, 



L f k * +> L 
E mf =-— (2t cos(fca) + 2t') 2 dk - — / (2t cos(fca) -2t')2 dk 

2tt y_ fe (+) 2?r J„ fe ( r ) 

4£ 9 8L , 2 

+— i 2 + — 1'\ (22) 

Notice that t < just inverts the cosine curves, translating the Brillouin zone considered in n/a, and t' < would 
just trade the roles of the two bands. Then, the relevant sector in the t, t' plane is t > and t' > 0. In this sector, 
the expression for the energy is 

— E mf = -8t sin(fc^T ) a) - ^'(tt - 2A4 _) a) + — i 2 + ^' 2 , f < t, (23) 

and 

a7r ^ / 47T 9 87T ,2 / 

—E mf = -Ant' + —t 2 + —t'\ t' > t. (24) 

The analysis of the above equations shows that, for J' < 8/ir 2 J, the MF configuration depicts two bands which 
coincide with those corresponding to two decoupled chains (t = J/tt, t' = 0), with Fermi momentum kp — kp ^ — 
k^ = 7r/(2a). Notice that the condition in eq. ( |l4] ) holds and the linearization procedure around this minimum is 
valid. 

On the contrary, for J' > \j , we find that the global energy minimum corresponds to the point t = 0, t' = J'/4 
where the condition in eq. ( |l4| ) does not hold (notice that there is still another local minimum while J' < 2J). The 
system in this configuration, which describes the strong coupling phase, presents a finite energy gap to spin excitations. 

In the following two sections we will explore the J' < 8 / it 2 J region. 



III. FLUCTUATIONS AND CONSTRAINTS: THE SU(4)i/U(l) COSET THEORY 

In this section we take the continuum limit of the MF Hamiltonian in eq. ( |l3| ) and then include fluctuations around 
MF and constraints in eq. (0). The outcome of this procedure is a perturbed SU(A)i/U(l) coset theory. 



A. Low-energy linearization in the thermodynamical continuum limit 

In the region we consider (J' < 8/ir 2 J) the mean- field dispersion relation consists of two coinciding bands of 
amplitude 2J/tt. Linearization of low-energy excitations can be done around kp = -k /(2a) in the usual way. The 
bandwidth will limit the validity of the resulting effective Field Theory to energies much smaller than J, independently 
of J'. 

Low-energy excitations in the thermodynamical continuum limit of the tight-binding model at half filling can be 
linearized in terms of Dirac fermions [p!l| . Fermionic position space operators for each band are readily written 
in terms of Dirac fermions (x) as 

= Va~(exp(-ikpx)i> ( +l(x) + exp(ik F x)^+l(x)), (25) 

<4;2 = Va{exp(-ik F x)ip { xl(x) + ex^(ik F x)^ L ' ) a (x)). (26) 

Here ^i^l and ^i 4 ^ stand for the right and left components of a Dirac spinor and so on. Dirac gamma matrices 
are taken as 70 = <ti, 71 = 02. Notice that there is a total of four Dirac fermion species; using the notation a =f, J, 
they are (+, f), (— , j), (+, J.), (— , 1), which will be respectively denoted ^(2:), with i = 1,2,3,4. Summation over 
repeated fermion species indices will be understood. 
We arrive then at the linearized MF Hamiltonian 

H mf = v F J dx^iix^d^iix), (27) 

where 

vp = 2tas'm(kpa) = 2Ja/ir (28) 

is the Fermi velocity, and $j = 5 , j7o- All four ^-fields have the same Fermi velocity, then the model, up to this point, 
possesses a manifest U (4) symmetry. 



5 



B. Fluctuations around mean field 



We include now the fluctuation fields B„\ B„\ B' n . As we look for the continuum limit of the Hamiltonian (|^), 
we will keep only relevant powers in a, as compared with eq. (p7|). 

In order to keep track of a orders, it is useful to make explicit the order a contribution of fermion bilinears by 
defining 

Ji) _ „-i„t(i) Ji) 

J2) _ 1 t(2) (2) ' / 2 q^ 
_(3) _ -1 t(l) (2) ' 

so that the leading order for each z is a . Notice that z^> and z^ still have to be expanded, as (n + l)a = x + a; the 
only relevant term in this expansion is that linear in a, containing first derivatives of "0-fields. Our notation will be 



= u,W + aw W (30) 



(notice that u„ = 0). 



The relevant expansions for the S-fields, taking into account that the MF value of B'q vanishes, are 

B^ =B + iaB A^ + aBoRiP + 0(a 2 ), {I = 1, 2) s 
B' n = iaA' n + aR' n + 0(a 2 ). { ' 



In particular, the terms quadratic in B must be expanded as 

= Bl + 2aB 2 R^ + 2 a 2 B 2 i?« 2 , (/ = i, 2 ) 



(32) 



B' r }B> n = a 2 A> 2 n + a 2 R>l. 

Using all of these, and making explicit the sums over I = 1, 2, the effective low energy Hamiltonian for (j^) is written 

as 

H eff = H {1) + H [2) + H {3) + 0{a 3 ) (33) 

where 

F« = J -NBl + aB ^ + 4 (1) ) + £ £ m 2 5o4 X) " + 



ra=l 



| £ a(aflbi#W + + Mgjtf' + 2a5 2 ii«*), (34) 



(2) _ 



-NB 2 + - 
2 2 



JV 



JV 



£as (4 2) + 4 (2) ) 



J 



N 



± £ aMoi?^^ 2 ) + W t< 2 )) + 25 2 i?( 2 ) + 2a£ 2 i? 



(2)' 



(35) 



V N 



2 

J' N 



^E a2 K(^i 3) +^ (3) ) + <] 06) 



2 

n=l 



The main things to notice here are: 

- there are irrelevant (divergent) constant terms. This is expected from the combination of Hubbard-Stratonovich 
and MF techniques. 
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- the terms without fluctuations in and provide the two decoupled chains MF results discussed in the 
previous section. 

- The AM and fields act as La grange multipliers; their total contribution to the effective action in the continuum 
limit reduces to a term 



dx (^ l {x) 11 ^. l {x))(A^{x)+A^{x)) + (^(x) 7 i(ai ® l)^ 3 (x))(A^ (x) - A^(x))) . (37) 



.Vf 



In the notation of eq. ( |37| ) the first matrix (<j\) refers to isospin indices (+), (— ), while the second one (1) refers to 
spin indices f , J. . 

- The presence of a quadratic term in the A 1 field, with proper sign, allows for a trivial Gaussian integration. The 
same can be done with the R fields. These of course brings back the original spin-spin rung interactions. In the 
present scheme their contribution includes quadratic terms in the c operators, that lead to a redefinition of the Fermi 
velocity vf —> w.f/2, and quartic perturbations that can be arranged as 



(38) 



The continuum form of these quartic perturbations in terms of Dirac fermions is lengthy. We will write them down 
below, after introducing a convenient notation. 

We notice that, for J' < 8/ir 2 J, our approach leads to a description of the system which is the same as the one 
obtained in perturbative treatments, in principle valid for J' -C J |4[JT^]. In particular, the first two terms in ( |38| ) give 
rise to the well known marginally irrelevant perturbation terms in the individual chains. However, our approach does 
not rely on any perturbative treatment of J' and in particular allows for the determination of the phase diagram of 
the system, i.e. it predicts a critical value of the ratio J'/ J which separates the two different regimes in the two- leg 
ladder. The situation is depicted in Fig. 2. 



J' 



I I I I A 



FIG. 2. Phase diagram of the spin ladder. Bold bonds correspond to non zero links in the MF approximation. 

Moreover, we show in the next section that the weak coupling structure unraveled in 0] arises naturally within our 
approach. 



C. Constraints 

In this section we express the constraints ([|) in terms of the linearized fermion fields and discuss how to implement 
them in the evaluation of the partition function for the spin ladder. 

In the continuum limit the constraint on the occupation number at each site (l),n separates in four parts, corre- 
sponding to oscillating and non-oscillating terms associated to each band. They read: 

^i'jo^i — constant (39) 
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which, implemented through a Lagrange multiplier ao, provides the time component of a gauge field implementing a 
diagonal (7(1) coset constraint, a M = (ao,A^ + A^); 

f i7o(oi ® = (40) 

which, implemented through a Lagrange multiplier bo, provides the time component of a second gauge field imple- 
menting the isospin U(i) coset constraint, — (bo,A^ — ^4(2)). These last two constraints 

= 0; (41) 



*i(<ri ® a)y*j =0 , (42) 
lead to marginally irrelevant quartic perturbation terms when implemented through 

S(0(x)) oc lim e -voHx)o(x)^ ^ 

t] — >oo 

just as in the case of decoupled chains. 

To conclude this section, we collect all the terms in the effective low energy Hamiltonian which finally reads 

H eff = ^ J dx &i(x) ((7i9 x - i-y^) 8 t] - i-y^ai ® 1);A) Vj(x) + AH eff , (44) 

where AH e ff includes quartic terms in fermionic fields, which arise from ( |38| ) and pi]) , ([i"2|) . 

Clearly, the unperturbed theory posses a 17(4) symmetry which is gauged by a diagonal (7(1) field (et p ) and an 
isospin (7(1) field (6 M ) which leads to the coset 

m si/®! (45) 



P(l)& 9 x[/(l) lM U"(l) iao ' 



where 517 (4) i stands for the level k = 1 WZW theory 

Before displaying the explicit expression for the perturbations it is worth discussing in more detail the coset structure 
of the quadratic part of the Hamiltonian. 



IV. 5!7(2)2 x Z 2 EMBEDDING, THE PERTURBATIONS IN A NEW LANGUAGE 

As it is known, the coset CFT 5f7(4)i can be alternatively described through the embedding ]lq ] 

5[7(4)x = SU(2) 2 x SU(2) 2 . (46) 

The conformal central charges of the two theories coincide and primary fields in the SU (4) i can be written in terms 
of primaries in the two SU(2)2 sectors. This will presently prove to be useful in the treatment of the perturbations. 
The different SU(2)2 sectors in this embedding are naturally identified in eq. ( f44| ) as the spin and isospin sectors, in 
virtue of the a\ ® 1 non-diagonal structure. Moreover, in this language the second Lagrange multiplier b^ gauges a 
(7(1) subgroup of the isospin 5(7(2)2 sector giving rise to 

SU{2)f m x SU(2) 2 sospm /U(l) lsospm = SU{2)f m x Z 2 . (47) 

The last factor has been identified in Q from the Z 2 structure of a two chain system. 

All of this is most easily shown in the bosonized version of the coset CFT. To this end we write fermion bilinears 
as Ull 

= M$g (48) 

where M is a renormalization constant and we have introduced bar indices in order to distinguish components 
transforming in the right and left fundamental representations of 5(7(4)i. The □ subindex indicates the fundamental 
representation in the standard Young tableaux notation. In identifying the two SU(2) 2 sectors we find useful to keep 
the original spin and isospin (band) indices, writing 
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3> 



aot,b(3 



(49) 



where we now use o, 6 for (+), (— ). 

This field $□ has scaling dimension 3/4 and its components can be written in terms of products of the components 
of the fields in the fundamental representations of the two SU {2)2 sectors as 



(50) 



where <p n and <fi' n are the primary fields in the fundamental (spin 1/2) representation of the two SU (2)2 isospin and 
spin sectors respectively. These fields have scaling dimension 3/8 so the product has the right dimension 3/4 and 
moreover, correlation functions of the fields on both sides coincide. 

The other primary field in the SU(A)i CFT is the one transforming in the antisymmetric (6 x 6) representation, 

which in the Young tableaux notation should read It is built up from the antisymmetric product of two fields in 

□ 

the fundamental representation 



□ 



.4 ($□$□) 



This field has scaling dimension 1 and can be mapped into SV (2)2 fields as 



(ft [(aiai),(a 2 ,a 2 )],[(bi/3i 
□ 



,(b2,02)] 



I {a 1 a 2 },{b 1 b 2 } ai a 2 ,0102 
!>□□ t t 



,'{aia2}{/3i/3 2 } 



(51) 



(52) 



where <pnn,4>'nu are the primary fields in the symmetric (spin 1) representations of the two SU (2)2 sectors, which 
have the correct scaling dimension 1. In eq. (|52| ) we have used the symbols {,} and [, ] to indicate respectively 
symmetrization and antisymmetrization of indices. 

We are now ready to analyze the different perturbation terms in AH e ff. First of all, contribution coming from 
intrachain couplings and constraints are known to be marginally irrelevant, just as in the case of decoupled chains 

The interchain perturbation terms in AH e ff (those arising from the last term in (|38|)) can be separated into 
two groups according to their scaling dimensions: there are terms which correspond to relevant operators (scaling 
dimension 1) which can be identified with certain linear combination of the components of the primary ( (5l| ) in the 
coset theory ([44] ) and current-current terms, which have scaling dimension 2 and are hence marginal. 

More precisely, for the relevant part we can write 



relevant perturbations = — A J dx ^Tr ^A$n^ + H.< 
where A oc J' and A is given by 



(53) 



(- 1 

1/2 





-1/2 

V 






-1/2 1/2 

1/2 -1/2 





\ 

-1/2 





1/2 

-1 J 



(54) 



(see the Appendix for details). Using the identifications described above and after some straightforward algebra we 
can readily identify the perturbation terms (p3[) in the embedding theory as 



(55) 



relevant perturbations = —A J dxTi ((f>nn + H.c.) + — J dxTv (^j}' n ai(/)^ai + H.< 



To analyze the effect of these perturbation terms it is convenient to reformulate the SU (2) 2 WZW sector in terms of 
three decoupled Majorana fcrmions, and in this new language it is easy to see that the first term gives a mass to all 
three Majorana fields |Tij] . The second one is simply the energy operator of the remaining Majorana sector 
Being all perturbations of dimension 1 we see that the gap opens linearly with the interchain coupling as predicted 
from the weak coupling limit ||,f|]. Note the different sign in the masses of the two sectors, also in agreement with 
the weak coupling analysis. 

As for the current-current terms, they correspond to marginal perturbations and can be written as 
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marginal perturbations = — — / dx ^— ^l((o"3 + ici) ® HOV'-k^lCO 7 ^ — io"2) ® — R ^ Lj 

J dx {i>\{(Tx ® l)Vi^(<Tl <8> + ^1(^1 ® &)iI>i.4r(<ti ® (56) 

The first two terms (first line) renormalize the Fermi velocity of the Majorana (isospin) sector, while the third and 
fourth (second line) correspond to marginal forward-scattering terms in the same sector. The fifth and seventh terms 
are effectively zero due to the constraints on the two corresponding U{1) currents and the sixth term correspond to 
the marginal forward-scattering terms in the spin sector. The very last one mixes spin and isospin sectors. This 
last contribution is nevertheless marginal, so it does not change the low energy physics which in the present case 
is dominated by the relevant perturbations already discussed. Its effect could be important in the analysis of e.g. 
zig-zag ladders where the relevant perturbations are wiped out, as we show in the next Section, and only marginal 
interactions play a role ]iT|-|22| . 

It can be easily shown that the marginal terms which are present on each separate chain, written in the present 
language correspond to the sixth and eight terms in the above expression. Due to the fact that these terms correspond 
to marginally irrelevant couplings and that they form a closed algebra, they will have no effect in the low energy 
dynamics whatsoever. After having observed that, one can see that the effective theory consists of two sectors which 
are decoupled from each other. 



V. OTHER STRUCTURES: CROSSED AND ZIG-ZAG LADDERS 



In this section we will extend our previous analysis to more general situations, which are not only of academic 
interest, but are relevant in the analysis of real materials. These more general situations arise when other (diagonal) 
couplings between spins in neighbouring chains are not negligible. The two structures that we analyze now are the 
so-called crossed ladders |2^,|2f|, in which couplings along the two diagonals are added, and zig-zag ladders in which 
only one diagonal coupling is added |l9|-|2^|. Another potential application of the present formalism would be the 
study of the interplay between interchain coupling and dimerization along the legs [p6|-]2§|] . 



i) Crossed ladders 



We consider a Heisenberg Hamiltonian given by 

N 2 N N 

H = E E JWJ&i + J' E + J *E (^■S'nli + sUi-SP) ■ (57) 

n=l 1=1 n— 1 n— 1 

where the last term corresponds to additional diagonal couplings. 

Following the same approach as for the normal ladder we introduce Hubbard-Stratonovich fields associated to each 
coupling and perform a three-parameter MF analysis proposing constant values for the intrachain couplings, the 
interchain (rung) coupling and the interchain (diagonal) couplings. We find two different regions in the parameter 
space (J/ J', J x /J'). It should be noted that this Hamiltonian is dual under the interchange J <-> J x , then it is enough 
to study the region J x < J. 



(a) If J/ J' > 7r 2 /8, the MF analysis yields the system in a weak coupling regime, and following all the same steps 
as before, we arrive at the same effective Field Theory with the noticeable change that the coefficient of the relevant 
perturbations is now shifted as J' — 2J X . As in the weak coupling analysis (2^j24|], one immediately sees that there 
is a line in which the relevant perturbations vanish. On this line one could expect a massless regime, as suggested by 
numerical studies [^3 25 1. However in a recent treatment of the resulting bosonized Hamiltonian it was shown that the 
current-current terms are marginally relevant and a gap opens |24}| . The same conclusion is attained in our resulting 
effective theory. Again, the new feature here is that we find the region of validity of the weak coupling effective Field 
Theory to go up to J' = 8/tt 2 J. 
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(b) If J I J' < 7r 2 /8, the system falls in a strong coupling regime in which the two dispersion bands are separated by 
a gap (oc J') and then a low-energy effective Field Theory description is not suitable here. 



j iy 
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XIX 



K 2 
8 



J/J' 



FIG. 3. Phase diagram of the crossed ladder. Bold bonds correspond again to non zero links in the MF approximation. 



ii) Zig-zag ladders 



The Hamiltonian is given by 



JY 



H — J^n "* 'Sn+1 



n=l 1=1 



N 
n=l 



(58) 



Introducing again Hubbard-Stratonovich fields associated to each coupling and performing a MF analysis with constant 
values for the intrachain and interchain couplings, we find a different situation: while we still find a regime, which 
now exists for J' < J, in which we re-obtain the standard weak coupling results, we find that the "strong coupling" 
regime, (J 7 > J), can still be described by an effective low energy Field Theory. 

More precisely, in the regime in which J' < J we find that all relevant perturbations cancel in a way similar to 
that found in the weak coupling limit 19-221. The effective low energy theory cor resp onds to the same coset theory, 
perturbed only by the operators appearing in the first, third and fourth lines in eq. J5q). The so-called parity breaking 
terms first studied in pjj ] appear in the pre sent approach from the next-to-leading order in the lattice-spacing a in 
the expansion of the modified version of ( |38|) . 

In the other regime, (J' > J), the bands at the MF minimum are given by 



,(+) 



a/2(1 + cos(fca)), 



J' 



(k) = — X2(l + cos(/fca)) 

7T 



(59) 



being no gap between them, and a Field Theory description is still possible. The difference is that the low energy 
effective theory should in this case be built up on only two fermion species, exhibiting SU(2)i symmetry. This should 
correspond to the description of a single chain plus next-nearest-neighbour interactions, which is the suitable picture 
for the regime where J' dominates. 

Once again, our method allows for the construction of an effective field Theory for the full range of couplings and in 
particular would allow to study the transition from the massless (c = 1) J = limit to the massive Kosterlitz-Thouless 
regime known to arise at 0.24 J' |2{|, which should, according to our analysis, extend to the limit J' — > 0. 

Since the main purpose of the present paper is to emphasize the potential applications of our approach, the analysis 
of these effective field theories will be addressed in a separate publication. 
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FIG. 4. Energy bands for the zig-zag ladder at J' > J MF minimum 




FIG. 5. Phase diagram of the zig-zag ladder. Bold bonds correspond again to non zero links in the MF approximation. 
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VI. CONCLUSIONS 



The approach developed in the present paper shows that the spectrum predicted from weak coupling approximation 
extends up to a finite value of J'/ J, our estimation of this critical value being ~ Beyond this value our MF 
analysis of Section II predicts a cross-over to the strong coupling regime, where the rungs of the ladder become 
disconnected among them. Fluctuations over this state will restore connectivity and the strong coupling approach 
of |50|,|l] would be the appropriate starting point in this parameter regime. As the classical potential analyzed in 
Section II has a double well structure in the intermediate region (J 7 ~ J' c ) we expect a smooth cross-over from weak to 
strong coupling regime. Experimental observation of this cross-over supposes the variation of the ratio of the exchange 
parameter. This could in principle be achieved by applying pressure in the perpendicular direction of the ladder axis. 

Though our approach starts from a MF analysis, fluctuations are taken into account to all orders. Besides, it allows 
for a classification of all the perturbations in the language of the embedding of the theory into SU (2)a x Z%. One 
interesting observation which arises is that only the Z2 Majorana Fermi velocity is renormalized to first order in J' 
by the interacions. 

The study of hole doped spin ladders is a natural extension of our approach. For this case the t — J model 
should be considered and the charge sector of the theory could be represented by a spinless boson (the slave boson 
representation). However the magnetic excitations will evolve from the triplet and the singlet found in this paper. 
The question of the hole pairing due to these excitations could therefore be addressed within our formalism. This will 
be reported elsewhere. 
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APPENDIX 



We write in this appendix the explicit form of some lengthy expressions appearing with compact notation in the 
main text. 

The relevant part of the third term in eq. (J3q) , appearing in eq. (H3) , reads in the continuum limit 



relevant perturbations = 
J'a 



(+) 

L,0 



— / dx (_tft(-)^t(-)^(-)^(-) , *t(-)^t(-)^(-)^(+) _ #tj-)*t(-U(+)^(-) , ^t(-)^t(-U(+)^ 

2 / UJ/ \ *R,a * R,0 w L.u* LJ3 T * R,a ^ R.J3 * L,q w L.fi w R,a * R,0 * L,a w L,0 T * R,a w R,/3 ^L,ct 
^R,a * R.p ^L,a^L.0 + * R.a W R,p ^L,a^L.f3 * R,a W R,0 ^L,a^L.0 + * R.a W R,p W L,a^ L,0 



by simple use of eqs. (|25|), (j26|), (|29|), (J30|) - 

The explicit form of eq. ( pl[ ) in terms of fermions, using eq. (fl8|), is 

<jd(ilo:i), (12,0:2)], [Cji0l),(Ja,#O] _ ( < j ) ( ] n"i).(Ji/3i) ( j)fe,a2),(j2,fc)\ 

where antisymmetrization affects bar and unbar pairs of indices separately. 

Using eqs. ( |60|) and (|6l[), expression ( |53| ) follows immediately. The base used for writing the matrix A in eq. ( |54| ) 
is the one made explicit with indices in the l.h.s. of eq. (p2|), ordered as [(+, |), (— , f)], [(+, T)> (+> !)]> [(+> T); ( — > I 

)],[(- T), (+,!)], [(-,!),(-, 4)], [(+,4), (-,!)]• 
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